Negative length orbits in normal-superconductor billiard systems 
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The Path-Length Spectra of mesoscopic systems including diffractive scatterers and connected to 
superconductor is studied theoretically. We show that the spectra differs fundamentally from that 
of normal systems due to the presence of Andreev reflection. It is shown that negative path-lengths 
should arise in the spectra as opposed to normal system. To highlight this effect we carried out both 
quantum mechanical and semiclassical calculations for the simplest possible diffractive scatterer. 
The most pronounced peaks in the Path-Length Spectra of the reflection amplitude are identified 
by the routes that the electron and/or hole travels. 
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In recent years, semiclassical methods have become 
a popular tool for describing devices operating in the 
mesoscopic regime. Advances in manufacturing and ma- 
terial design have made possible the creation of clean 
mesoscopic devices, whose properties depend on the mi- 
croscopic details of individual samples. For example, in 
recent experiments jl],^| involving semiconductor micro- 
junctions, both the quantum coherence length and the 
mean free path of elastic collisions are large compared to 
the size of the junction. In such devices electrons can be 
described as a two-dimensional ideal Fermi gas of non- 
interacting particles. The conductance of such junctions 
has been measured and found to oscillate strongly as the 
Fermi energy is varied. 

Semiclassical methods have proved to be very effec- 
tive for understanding conductance fluctuations in nor- 
mal microjunctions. On the one hand, methods based 
on random matrix theory successfully predict sta- 

tistical properties of transport properties. On the other, 
short-wavelength semiclassical descriptions able to ex- 
plore geometry-induced interference effects in weakly dis- 
ordered or clean mesoscopic devices. In particular, it has 
been shown that the Path-Length Spectra (PLS), defined 
as the power spectrum of the reflection (transmission) 
amplitudes with respect to the Fermi wavelength 
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possesses peaks at lengths corresponding to classical tra- 
jectories of electrons starting and ending at the external 
contacts ]§-H]. Here r mn (kF) (t m n(kp)) is the reflection 
(transmission) amplitude at the Fermi wavenumber k F 
for scattering from mode n of the entrance lead to mode 
m of the entrance (exit) lead in a two probe conductance 
measurement. 

Mesoscopic devices connected to a superconductor 
present a new challenge for semiclassics. In such 



normal-superconductor (NS) systems Andreev reflection 
]i~3|,p| p4p^ | plays an important role, whereby electrons at 
the Fermi energy in the normal metal are retro-reflected 
as holes at the NS interface. Such a process might be 
expected to dramatically affect the PLS, but to-date, no 
investigations of the PLS in NS systems have been carried 
out. 

In this Letter we present the first such investigation, 
by examining a mesoscopic device connected to a super- 
conductor and show that the PLS of NS systems differ 
from those of normal systems in a fundamental way. In 
particular for NS systems containing diffractive scatters, 
negative path-lengths can arise in the PLS which are ab- 
sent from the corresponding normal systems. 

Before turning to the NS system we shortly discuss the 
role of diffraction in the PLS of normal systems, following 
Rcf. [fl6|| . In the semiclassical approximation, particles 
hitting a diffractive scatterer may scatter in any direc- 
tion since the classical dynamics is not uniquely defined 
jl7],[l8]]. In Ref. |ll| a two-dimensional wave guide with 
a small point-like diffractive scatterer has been analyzed 
]l9[ (see Fig la.). It has been shown that the PLS of the 
reflection amplitude r nm has peaks at path-lengths corre- 
sponding to classical trajectories starting and returning 
to the entrance of the lead, either diffracted once or sev- 
eral times by the scatterer. At large Fermi wavelengths 
a typical trajectory with multiple bounces is shown in 
Fig. la. Such trajectories consist of two parts: segments 
(along the z axis) connecting the lead and the scatterer 
with total length 2 x zq, and multiple diffraction trajec- 
tories starting and ending on the scatterer. A multiple 
diffraction trajectory can be decomposed into loops that 
start and end on the scatterer, making bounces on the 
walls of the waveguide. The possible lengths of the loops 
are 
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where r = 0,1,2,... is the repetition number and W is 
the width of the waveguide . 

In Figure || the PLS calculated quantum mechanically 
and semiclassically using the method developed in Ref. 
p6[ are shown. One can see that the agreement between 
the quantum and the semiclassical calculation is excel- 
lent and the peaks are located at lengths which are lin- 
ear combinations of (Q) plus 2zo- It is obvious that the 
PLS has peaks only for positive lengths L. The ampli- 
tude of the peaks is decreased by multiple diffraction and 
therefore the most pronounced peaks correspond to paths 
diffracted only once on the scatterer. 

Now consider the effect of replacing one of the exit 
leads by a superconductor. In this case a new contri- 
bution to the reflection amplitude r mn has to be taken 
into account, namely the one coming from Andreev re- 
flection at the NS interface. By solving the Bogoliubov- 
de Gennes equation for a NS interface it is possible to 
show that, at the Fermi energy, electron- like excitations 
impinging onto the superconducting interface are coher- 
ently retro-reflected as hole-like excitations. In a semi- 
classical description, the classical action associated with 
the path connecting points q' and q" is given by 



S{q',q") = / p(q)dq 



(3) 



where p(q) is the momentum of the electron or the hole 
along the path (note that if the path touches the super- 
conductor at least once, it will contain both electron and 
hole parts). In particular the action associated with an 
Andreev reflection process in which an electron starting 
at point q' returns to the same point as a hole can be 
written as 



S(q', q') = k F L(q', q*) - k F L(q*,q') = 0, 



(4) 



where L(q', q*) denotes the length of the path of the elec- 
tron until it hits the superconductor at q* and L(q*,q') 
is the path length of the hole from q* to q' . The minus 
sign in the second term is due to the fact that the direc- 
tions of the momentum and the velocity of the hole are 
opposite. The total action of an electron-hole (e-h) tra- 
jectory returning to its starting point is always zero, since 
the hole retraces the path of the electron. As a conse- 
quence of this result, the PLS of the reflection amplitude 
for electron-electron (e-e) has peaks at positive lengths 
L while for e-h it has a pronounced peak at L = 0. 

We now consider the case where diffractive scatter- 
ing is possible in the system (Fig. lb). The new fea- 
ture is that the trajectory of electrons or holes hitting 
a diffractive point is not uniquely defined. Therefore, a 



hole retracing an electron, which scattered on a diffrac- 
tive center, will not necessarily retrace the trajectory of 
the electron beyond the diffractive center. The hole may 
leave the diffraction center at a different angle than that 
of the incident electron. This effect has already been 
pointed out by Beenakker pTj| in connection with normal- 
metal-superconductor junction containing a point con- 
tact. Consequently, in the presence of diffraction, compli- 
cated trajectories consisting of several electron and hole 
segments may arise. The classical action for this case can 
be written as a sum of actions of the segments 



S = ^2k P (±Li), 



(5) 



where Li is the length of the segment i and the + or — 
correspond to cases when electron or hole traverses the 
segment i, respectively. Unlike in the diffractionless case, 
the sum of positive and negative terms in Eq. (||) is not 
necessarily zero. Moreover, the total length of the hole 
segments may exceed those of the electrons. Thus, in the 
PLS of the e-e reflection amplitude, peaks at negative 
lengths may appear, which are completely absent from 
the PLS of the corresponding normal systems. 

To observe negative lengths in the PLS we note that 
when the exit lead is replaced by a superconductor (|] 
the e-e reflection amplitude can be expressed in terms of 
the transmission and reflection amplitudes of the corre- 
sponding normal system (see Eq. 266a in Ref. ||). For 
the system of Fig. lb, the latter amplitudes have been de- 
rived in Ref. fl6|| and, after lengthy but straightforward 
algebra, we find the following expression for the n, m ma- 
trix element of the e-e reflection amplitude at the Fermi 
energy: 
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2r nm (k F ) 



1 + \V\ 2 [ImG (x ,z \x , z )Y 



(6) 



where r nm (k F ) is the matrix of reflection amplitudes 
of the normal system |TJj for the entrance lead, 
Go(x, z\x' , z') is the Green's function of the empty waveg- 
uide. T> = 1 — AGo(a;o, zoNch zq) can be regarded 

as the diffraction constant of the scatterer, where A is the 
renormalized strength of the scatterer |2l| . 

On the upper part of Fig. || the PLS of the reflection 
amplitude of (^|) (the exact quantum mechanical result) 
is plotted as a function of the path length. Peaks at 
negative path lengths are clearly visible here as opposed 
to normal system shown in Fig. |^. 

The semiclassical approximation of Eq. (Jg|) can be ob- 
tained from the semiclassical form of the Green's function 
which is given by 



G (xo,zo\x ,z ) = ^2 



Q ikFl r — i3-K /4 



(7) 



where the summation is over all possible loops with 
lengths l T given in (0) and n r is the number of bounces 
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on the walls of the waveguide. One can see that 
ImGo(xo, zq\xq, zq) will contain the complex conjugate 
of Eq. (M) and is therefore a sum with terms proportional 
to e ±lk ^ lr . Here, terms with e ~ lkFlr correspond to loops 
traversed by holes. The c-e reflection amplitude given in 
(Q) can be rewritten as a multiple diffraction scries: 

C^f) = 2r nm (k F ) jl - \V\ 2 [Im G (x , z \x , z )f 
+ \V\ 4 [ImGo(x ,z Q \x Q ,z Q )] 4 ---y (8) 

Expressing powers of Im Go with Eq. (Q) and its complex 
conjugate s^^ikp) can be written as an oscillating sum 
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where Sj is of form Eq. (||). The amplitudes Aj can 
be determined exactly by using the above formulas. In 
the lower part of Fig. || the PLS of the reflection ampli- 
tude s^(kF) is plotted using Eq. (||). Again one can 
see a very good agreement between the quantum and 
semiclassical calculations. The most pronounced peaks 
with negative lengths (see Fig. |^) come from the family 
L = 2zq — l r , where l r is given in Eq. (^), namely L = 
-0.4, -1.0,-1.6, -2.4, -3.0, -3.6, -4.4, -5, ■ • • . These 
lengths can be associated with the following routes: First 
the electron hits the superconductor then the retro- 
reflected hole diffracts on the scatterer and makes a loop 
with repetition number r as described before Eq. (||). 
Next, the hole diffracts off the scatterer then hits the su- 
perconductor on which it converts back to an electron. 
Finally, the electron goes back to the entrance lead. 

The above results represent the first theoretical study 
of the path-length spectra of a normal superconductor 
mesoscopic system with diffractive scattering. We have 
demonstrated that the PLS of such systems differs funda- 
mentally from that of normal systems, due to the appear- 
ance of peaks at negative lengths. To highlight this effect, 
we have analyzed the simplest possible diffractive scat- 
terer. Since the appearance of negative lengths in PLS is 
the direct consequence of the presence of diffractive scat- 
terers it is desirable to study other types of diffractive 
centers such as corners. There is a growing interest in 
studying the role of diffractive scatterers in normal meso- 
scopic systems p2] |, therefore the extension to normal 
superconductor systems may become a new playground 
both in the semiclassical theory of scattering processes 
and level statistics of these systems. For the future it 
will also be of interest to examine the amplitudes of the 
negative-length peaks in more complex geometries, such 
as those of |15j, to examine the effect of a tunnel junc- 
tion at the NS interface and the role of order parameter 
symmetry. 
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FIG. 1. The waveguide with point-like diffractive scatterer 
(a) opened and (b) a superconductor attached to the exit 
lead (left side of the waveguide). The origin of the coordinate 
system (z, x) is at the left bottom corner of the waveguide. 
In Fig. la one possible trajectory is shown. For clarity the 
forward and backward paths are shifted. 
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FIG. 3. The PLS of the reflection amplitude s 11 " e (fc F ) (see 
Eq. (|^)) is obtained from quantum and semiclassical calcula- 
tions. For parameters see Fig. 2. 
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FIG. 2. The PLS of the reflection amplitude defined by 
Eq. ([j]) is obtained from quantum mechanical and semiclas- 
sical calculations. The width of the waveguide, W = 1.0. 
The strength of the Dirac delta potential, A = 10.0, and its 
position is located at {zo,xo) = (0.5,0.7). 




FIG. 4. The PLS of the reflection amplitude rn(k F ) 
of the normal system (lower part) and s^ e (kF) given 
in Eq. (^) (upper part) are plotted. The most pro- 
nounced peaks with negative lengths come from the fam- 
ily L = 2zo — l T , where l r is given in Eq. (Q), namely 
L = -0.4, -1.0, -1.6, -2.4, -3.0, -3.6, -4.4, -5, • ■ ■ . For pa- 
rameters see Fig. 2. 
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